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What is claimed is: 

1. A method for simultaneously determining respective scale factoraror 
alignment angles of sensitive axes in a multi-axis accelerometer device^or 
measuring acceleration, comprising the steps of: 

a) mounting a multi-axis accelerometer device on a turntable in a first 
orientation, the turntable having a tilt angle with respect to a/vertical axis 
defined by a local gravity vector; 

b) spinning a multi-axis accelerometer device aradnd an axis of rotation 
at an angular velocity using the turn table such that/oie multi-axis 
accelerometer device experiences a time varying component of the local 
gravity vector; 

c) receiving respective outputs of the^multiple axis as the multi-axis 
accelerometer device experiences the tinje varying component of the local 
gravity vector; 

d) repeating steps (a), (b) anflf (c) with the multi-axis accelerometer 
device mounted in a second orientation; and, 

e) repeating steps (a), (tO and (c) with the multi-axis accelerometer 
device mounted in a third orientation; and, 

f) determining respective scale factors or alignment angles of the 
multiple axes of the accelerometer device by combining the respective received 
outputs of the accelerometer device with predicted outputs of an ideal 
accelerometer, the$redicted outputs based on the tilt angle of the turntable, 
the angular velocity of the ideal accelerometer, and the local gravity vector. 
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5 2. The method of Claim 1 wherein the angular velocity is constant 

during the receiving. 
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3. The method of Claim 1 wherein the multiple-axis accelerom^ter 
device is oriented in three orientations while recording data. 

4. The method of Claim 1 wherein the time varying ownponents of the 
local gravity vector are equal to g*sin(9)*cos(<j)(t)) and g*&in(0)*sin(<|)(t)), where 
0 is the tilt angle, g is the acceleration due to gravity, and <J) is an angle 
subtended at the axis of rotation by the accelerome^er and the component of 
gravity in the plane of rotation of the accelerometer. 

5. The method of Claim 1 further including the step of filtering the 
outputs of the multiple axis using respegttve low pass filters. 




6. The method of Claim 5 further including the step of sampling the low 
pass filtered outputs of the multiple axis using respective analog to digital 
converters. 



7. The method of Claim 6 further including the step of receiving the 
25 sampled outputs of the imiltiple axis and combining the sampled received 
outputs of the multiple/axis with one or more predicted outputs to determine 
the scale factors of tHe sensitive axes. 



8. The method of Claim 6 further including the step of receiving the 
30 sampled outputs of the multiple axis and combining the sampled received 
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5 outputs of the multiple axis with one or more predicted outputs to determine 
the alignment angles of the sensitive axes. 

9. The method of Claim 1 farther including the steps of: 
taking respective Fourier transforms of the received outgjdts of the 
multiple axis; 

taking the Fourier transform of the predicted outpuj^ of an ideal 
accelerometer; and 

combining the respective Fourier transforms of the received outputs and 
the predicted output to determine the scale factors^r alignment angles of the 
0 15 multiple axis of the multi-axis accelerometer device. 




10. A system for simultaneously determining respective scale factors or 

4 alignment angles of a multi-axis accelerometer device for measuring 

3 ' 

; p acceleration, comprising: 

i ii 

20 a turn table mechanism configured to mount an accelerometer device 

having multiple axis for calibration, the turntable having a tilt angle with 
respect to a vertical axis defined by a local gravity vector, the turntable 
configured to spin the accepter device around an asis of rotation at an 
angular velocity such thatfthe accelerometer device experiences time varying 
25 components of the local gravity vector; and 

a processor system coupled to receive respective outputs of the multiple 
sensitive axes of the^ccelerometer device, the processor system configured to 
record the outputsZf the accelerometer device as the device experiences the 
time varying components of the local gravity vector and to determine 
30 respective scaWfactors or alignment angles of the multiple axis of the 

accelerometer device by combining the logged outputs of the accelerometer 
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device with a predicted output of an ideal accelerometer, the predicted 01 
based on the tilt angle of the turntable, the angular velocity of the ideal 
accelerometer and the local gravity vector. 

11. The system of Claim 10 wherein the turntable is cqjffigured to 
maintain a constant angular velocity during the recording. 



12. The system of Claim 10 wherein the time/varying components of 
the local gravity vector are equal to g*sin(0)*cos(4>it)) and g*sin(9)*sin((})(t)), 
where 9 is the tilt angle, g is the acceleration dp^to gravity, and $ is an angle 
subtended at the axis of rotation by the accelerometer and the component of 
gravity in the plane of rotation of the accelerometer device. 



13. The system of Claim lO^irther including a low pass filter for 
filtering the outputs of the accelerometer device. 



14. The system of Claim 13 further including an analog to digital 
converter for sampling tlje low pass filtered outputs of the accelerometer 
device. 

15. The system of Claim 14, wherein the processor system is further 
configured to determine the scale factors or alignment angles of the 
acceleromete^jdevice by recording the sampled outputs of the accelerometer 
device, and by combining the sampled, recorded outputs of the accelerometer 
device wit]/ the predicted output of an ideal accelerometer. 
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16. The system of Claim 15 wherein the processor system is furth^ 
configured to determine the scale factors and/or alignment angles of thf 
accelerometer device by: 

taking respective Fourier transforms of the recorded outj^ts of the 
multiple sensitive axes; 

taking the Fourier transform of the predicted output of an ideal 
accelerometer; and 

combining the respective Fourier transforms of the recorded outputs 
and the predicted output to determine the scale factors or alignment angles of 
the multiple sensitive axes of the multi-axis accelerometer device. 

17. A method for simultaneously determining respective scale factors or 
alignment angles of sensitive axes in a multi-axis accelerometer device for 
measuring acceleration, comprising the steps of: 

a) mounting a multi-axis accelerometer device on a turntable in a first 

injfa til 

defined by a local gravity vectofr; 

b) spinning a multi-axis accelerometer device around an axis of rotation 
at an angular velocity using the turn table such that the multi-axis 
accelerometer device experiences a time varying component of the local 



orientation, the turntable having^ tilt angle with respect to a vertical axis 



gravity vector; 



/ 



c) receiving respective outputs of the multiple axis as the multi-axis 

/ 

accelerometer device experiences the time varying component of the local 
gravity vector; / 

d) det^mining respective scale factors or alignment angles of the 
multiple axe^ of the accelerometer device by combining the respective received 
outputs of tihe accelerometer device with predicted outputs of an ideal 
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accelerometer, the predicted ojaxputs based on the tilt angle of the turntable, 
the angular velocity of th^raeal accelerometer, and the local gravity vector. 

18. The method of Claim 17 further including the step of repeating steps 
(a), (b) and (c) witfn the multi-axis accelerometer device mounted in a second 
orientation. / 

19/f The method of Claim 18 further including the step of repeating steps 
(a), (b)yand (c) with the multi-axis accelerometer device mounted in a third 
orientation. 
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Appendix A: Computing Scale Factors and Alignment Vectors 

This appendix shows a solution method for the system of equations derived in 
the main body of the patent, to obtain the scale factors and the alignment 
vectors of the sensitive axes of the multi-axis accelerometer device. 



In Orientation 1: 

P A ,i= K*a A *g*sin(e)*exp(i*<t>,)*( A,- i*Ap (1) 

p b,i= K*a B *g*sin(e)*exp(i*<)>,)*( B x - i*B y ) (2) 

P c .i= K*a c *g*sin(6)*exp(i*(t) 1 )*( C x - i*C y ) (3) 

In Orientation 2: 

P A2 = K*cc A *g*sin(e)*exp(i*(|> 2 )*( A, + i*AJ (4) 

P B 2 =. K*a B *g*sin(e)*exp(i*()) 2 )*( B x + i*B 2 ) (5) 

P c 2 = K*a c *g*sin(6)*exp(i*<t> 2 )*( C x + i*C z ) (6) 

In Orientation 3: 

P A3 = K*a A *g*sin(0)*exp(i*(|)3) :i: ( A, + i*Ay) (7) 

P B .3= K*a B *g*sm(9)*exp(i*(|>3)*( B z + i*B y ) (8) 

P C 3 = K*a c *g*sin(e)*exp(i*^ 3 )*( C z + i*C y ) (9) 

Ideal Accelerometer with sensitive axis parallel to the plane of rotation: 

P_nominal = K*a_nominal*g_nominal*sin(G_measured) (10) 



Additionally, the peak DFT values are known from the recorded data and 
generated ideal data. 
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In the initial stage, the value of k can be determined from equation (10), by 
substituting the known values for Pjnominal, a_nominal, gjnominal and 
6_measured. 

In the next stage, the absolute value of a A is calculated. This is the scale factor 
of accelerometer A. By squaring (1), the following equation is obtained: 

I P A ,i 1 2 = I k 1 2 *(a A *g*sin(6)) 2 * | exptf*^) | 2 * | (A,- i*Ay) | 2 (11) 

•••|Pa,iI 2 = |K| 2 *(a A *g*sin(e)) 2 *{(A x ) 2 +(A y ) 2 } (12) 
Similarly, by squaring (4): 

|P A , 2 | 2 = iKl^Vsi^e)) 2 *^)^) 2 } (13) 
Also, by squaring (7): 

I P A>3 1 2 = |K| 2 *(a A *g*sin(6)) 2 *{(A Z ) 2 +(Ay) 2 } (14) 

Adding (12), (13) and (14) gives 

|P a , 1 | 2 +|Pa, 2 | 2 +|Pa.3| 2 = 

| k 1 2 *(a A *g*sin(6)) 2 *{(AJ 2 +(Ay) 2 +(A X ) 2 +(A z ) 2 + (A,) 2 +(Ay) 2 } (15) 

However, since A,, Ay, and Aj form an alignment vector [A^ Ay AJ, which is a 
unit vector, the following identity holds: 

[A, Ay AJ . [A^ Ay A z ] = 1 (16) 
This yields the identity 

(A x ) 2 +(A y ) 2 +(A z ) 2 =l (17) 

Substituting (17) into (15) gives 

I Pa,i 1 2 + I Pa* 1 2 + I Pa.3 1 2 = I k 1 2 *(a A *g*sin(6)) 2 *{2} (18) 
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5 The value ot A can be obtained by substituting the values for , P A 2 , P^ 3 and 
9, as well as the value for k from solving (10), into equation (18) and solving the 
equation (18) with the knowledge that a A is positive. 

The scale factors of accelerometer B and accelerometer C, which are a B and a c 
10 respectively, can be obtained in a similar manner. Thereby, the scale factors of 
the sensitive axes of the multi-axis accelerometer device are obtained. 



In the final stage the alignment vectors are computed. 
«3 Dividing (1) by (2) and rearranging gives: 



P A ,i* a B .*( B x - i*B y ) - P B>1 *a A *( A,- i*Ap = 0 (19) 
Similarly dividing (2) by (3) and rearranging gives: 

P B /a c *(C x -i*C y )- P Cfl *a B *(B x -i*B y ) = 0. (20) 

Dividing (4) by (5) and rearranging gives: 
20 P A2 *cx B *( B x + i*B z ) - P B /a A *( A, + i*A,) = 0 (21) 

Similarly dividing (5) by (6) and rearranging gives: 

P B2 *a c *( C x + i*C z ) - P c /a*( B x + i*B z ) = 0 (22) 



Dividing (7) by (8) and rearranging gives: 
25 P^Ob *( B 2 + i*B y ) - P B , 3 *a A *( A, + i*Ay) = 0 (23) 

Similarly dividing (5) by (6) and rearranging gives: 

Pb,3*cc c *( C x + i*C 2 ) - P c /a B *( B x + i*B z ) = 0 (24) 

Additionally, [A^A^AJ, [B x B y B J and [C^C y CJ are unit vectors. Rearranging 
30 (17) gives the following equation: 

(A x ) 2 +(A y ) 2 +(A z ) 2 -l = 0 (25) 
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Similarly, the following equations can be derived: 
(B x ) 2 +(B y ) 2 +(B 2 ) 2 -l = 0 
(C x ) 2 +(C y ) 2 +(C z ) 2 -l = 0 



(26) 



(27) 



The set of equations, (19) through (27), can be solved for the alignment vectors, 
[A, Ay AJ, [B x B y BJ and [C^C y C z ], using a multidimensional Newton- 
Raphson solution method. The nominal values of the alignment vectors can be 
used as initial values for the Newton-Raphson method. 

Thus, the values of the scale factors and the alignment vectors of the sensitive 
axes of the multi-axis accelerometer device can be obtained. 
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